Kinetics of Coalescence, Annihilation, and the q-State Potts Model in One Dimension 
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Abstract: The kinetics of the q-state Potts model in the zero temperature limit in one dimension 
is analyzed exactly through a generalization of the method of empty intervals, previously used for 
the analysis of diffusion-limited coalescence, A + A — > A. In this new approach, the q-state Potts 
model, coalescence, and annihilation (A + A — > 0) all satisfy the same diffusion equation, and differ 
only in the imposed initial condition. 
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The kinetics of the q-state Potts model has drawn 
much recent attention, in particular regarding its persis- 
tence probabilities [1-10]. An understanding has emerged 
that, as far as persistence is concerned, the kinetics of the 
q-state Potts model in the zero temperature limit, with 
Glauber dynamics and in one dimension, is equivalent 
to that of the Ising model, when a fraction 1/g of the 
spins are initially up [6]. It is also understood that if 
one regards the boundaries between domains as particles 
(kinks), then q = 2 (the Ising model) corresponds to the 
diffusion- limited annihilation process A + A — > 0, while 
the limit q — > oo corresponds to coalescence, A + A^ A. 
Here we extend the formalism of empty intervals, pre- 
viously used for the exact analysis of A + A — > A [11], 
to annihilation, by considering the probability that an 
interval contains an even number of particles, and we ap- 
ply this new approach also to the g-state Potts model. 
This results in one simple diffusion equation describing 
all three models (annihilation, Potts model, and coales- 
cence): they differ only in the initial conditions they im- 
pose. 

Consider diffusion-limited annihilation, A + A — > 0, in 
a one-dimensional lattice of lattice spacing a. The par- 
ticles hop randomly to the nearest site to their right or 
left, at equal rate T, and annihilate immediately upon en- 
counter. Let G n (t) be the probability that an arbitrary 
segment of n consecutive sites contains an even number 
of particles, at time t. (We assume that the system is 
infinite and homogeneous.) A site can be either empty 
or occupied by a single particle, so the particle concen- 
tration is 



c(t) = 



l-Gx{t) 



(1) 



Because the reaction A + A — > conserves the num- 
ber of particles modulo 2, the only way that the parity 
within an n-segment is changed is when particles at the 
edge of the segment hop out, or when particles outside of 
the segment hop in. To describe these events, we require 
F n (t) — the probability that an n-segment containing an 



even number of particles is followed by the presence of a 
particle at the (n + l)-th site. This can be expressed in 
terms of the G n (Fig. 1): 



2F n (t) = (1 - Gi) + (G n - G n+ i) 



(2a) 



Likewise, H n (t) — the probability that an n-segment 
containing an odd number of particles is followed by a 
particle at the (n + l)-th site — is 



2H n (t) = (l-Gx)-{G n -G n+l ) 



(2b) 



G„ +1 {" 



Fig. 1: Computation of F n (t). Empty (solid) rectan- 
gles symbolize n-segments with an even (odd) number 
of particles. Empty (solid) circles represent empty (oc- 
cupied) sites. The four events listed may be grouped in 
two different ways, yielding eq. (2a). One can also see 
that the probabilities of the first two events add up to 
F n + H n = 1 - G u yielding eq. (2b). 

The evolution equation for G n is now readily obtained: 



2r(f 1 rl _i — H n ~i + H n — F n ) 



where the first term on the r.h.s expresses the event that 
a particle at site n jumps out of the segment, leaving 
an even number of particles in the remaining n — 1 sites 
(and hence in the n-segment); the second term denotes 
the same case, but when there are initially an even num- 
ber of particles in the n-segment (that is, an odd number 
in the (n — l)-segment); and the last two terms pertain 
to a particle just outside of the n-segment, at site n + 1, 
jumping in. The factor of 2 accounts for the fact that 
each event may occur at the left or right edge of the 
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segment, with equal probability. Using eqs. (2), this be- 
comes 

%-G n (t) = 2r(G„_! - 2G„ + G n+1 ) . (3) 
at 

The case of n = 1 requires a special equation, since Go 
is undefined. Taking into account all the ways Gi might 
change, one finds 

^G 1 (i) = 2r(l-2G 1 + G 2 ) . 

Thus, eq. (3) may be understood to be valid for all n > 1, 
provided that one uses the boundary condition 

G (t) = 1 . (4a) 

Additionally, since the G„ are probabilities, we have 

< G„(t) < 1 . (46) 

We now turn our attention to the q-state Potts model 
in one dimension, at the zero temperature limit and with 
Glauber dynamics. (This is also known as the voters 
model, with q opinions.) Each site can assume one of 
q states, si, S2,...,s g . The model possesses q equiva- 
lent ground states, where all the sites are in the same 
state. If the system is started away from equilibrium, 
then neighboring sites of differing states react according 
to the scheme: 

J SiSi rate Y , , . 

{Sj Sj rater. [t>) 

Hence the interface between domains of different states 
may be regarded as particles (kinks) that perform regular 
random walks. Let the kinks between a domain of state 
Si on the left, and a domain of state Sj on the right, be 
denoted by Aij. The dynamics (5) implies the following 
immediate reactions between colliding kinks: 

Aij + A jk -> A ik , (6a) 

A^ + Aji -> . (6b) 

Notice that the inner indices (j , in the cases above) must 
be the same for kinks to meet. 

Focusing now on the kinks, and assuming an infinite 
homogeneous system, let G n (t) be the probability that 
an arbitrary n-segment contains kinks such that their to- 
tal number of indices of each of the q kinds is even, at 
time t. Because the reactions (6) conserve the parity of 
the numbers of indices, reactions within an n-segment 
do not affect G n . Thus, G„ changes only through kinks 
that hop into or out of the segment, exactly as in the 
case of the annihilation process. (Note that a kink must 
have different indices, so the parity is always affected by 
such transitions.) It follows that the G„ obey eq. (3), 
with the boundary conditions (4). Finally, if one merely 



counts the kinks, ignoring their indices, then their con- 
centration is given once again by eq. (1). 

For q = 2 the process (6a) never takes place (it re- 
quires at least three states) and the problem reduces to 
that of the annihilation process, A+A — ► 0. On the other 
hand, the probability that the sites to the left and right 
of the two kinks are in the same state Si tends to zero as 
q — ► oo. In that case process (6b) never takes place and 
the model is then analogous to coalescence, A + A — > A. 
For 2 < q < oo kinks may either coalesce or annihilate, 
and we have a mixture of the two processes. Annihilation, 
coalescence, and the q-state Potts model are all described 
by the same (discrete) diffusion equation, eq. (3), and the 
same boundary conditions, eqs. (4). 

The difference between the three models is in their im- 
plied initial condition. Assume, for example, initial con- 
ditions where each site is randomly in one of its q possible 
states. The probability that an n-segment contains even 
numbers of indices equals the probability that the first 
and last site in the segment are in the same state. For 
the case of initially random states this yields 

G„(0) = 1 . (7) 

While the initial state above is natural, it is by no means 
unique. Consider, for example, the initial state of Fig. 2a, 
for which 

[ ^ < n < to , 
G„(0) = < l ™ kqm < n < (fcq + l)m, k = 1, 2, . . . 
1 otherwise. 

(8) 

This corresponds to kinks that are initially evenly spaced, 
to sites apart from each other. Notice, however, that the 
initial states in between the kinks could be chosen differ- 
ently (Fig. 2b), thus dictating a different form of G„(0) 
and different kinetics! 

3 I 1 I 2 I 3 ' 1 I 2 I 3 
(a) 

1 I 3 I 2 I 3 I 2 I i I 3 
(b) 

Fig. 2: Possible initial configurations for the q-state 
Potts model. Here the system consists of alternating do- 
mains of equal length to, which corresponds to equally 
spaced kinks (the vertical bars), m sites apart from each 
other. In (a) the states of the domains alternate peri- 
odically: 1, 2, . . . , q, 1, 2, . . . while in (b) the states are 
random (but do not repeat). In this illustration, q = 3. 

Eq. (3) may be solved by standard techniques, for ex- 
ample, by Laplace-transforming with respect to time, fit- 
ting an exponential solution to the resulting difference 
equation, and finally inverting the Laplace-transformed 
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solution [12]. With the boundary conditions (4) and the 
natural initial condition (7), one obtains: 

G n {r) = - + q — r n e-^I n (2r')^, 

q q Jo 2t 

where r = 2Tt and /„(•) is the modified Bessel function 
of order n [13]. In particular, the probability that a site 
contains a kink, is 



persistence, but the new formalism can also handle cor- 
related initial conditions, where the three models might 
differ in nontrivial ways. Moreover, the new approach 
can be extended to inhomogcncous situations, along the 
same lines as the method of empty intervals (IPDF) [14]. 

We thank L. Glasser and D. Kessler for edifying dis- 
cussions, and we gratefully acknowledge NSF support of 
this work (PHY-9820569). 



1-Gi(r) = 2_l[] - 



e--7 l( 2r')^). (9) 



Thus, with natural initial conditions, the density of kinks 
in the g-state Potts model has a simple dependence on 
q: it is exactly equal to the density of particles in the 
coalescence model, A + A — > A, times (q — l)/q. 

For completeness, we note that for low initial concen- 
trations of kinks, passage to the continuum limit is jus- 
tified. Putting x = na, G n (t) — > G(x,t), and T = D/a 2 
in eq. (3), and taking the limit a — > 0, one obtains the 
diffusion equation 



(10) 



The boundary condition (4a) becomes G(0,i) = 1, and 
the concentration is obtained from the continuum limit 
of eq. (1); c(i) = — dG/dx\ x= o- Consider, for example, an 
initial concentration cq of randomly placed kinks, where 
the states of the domains in between are random (but do 
not repeat). Following simple combinatorial arguments, 
it can be shown that this corresponds to the initial con- 
dition 

G(x,0) = - + S—l e -^ C0X . 

q q 

Solving for the concentration, we obtain 



c(t) — coe z erfc(z) , z 



q-1 



which has the expected long-time asymptotic behavior 
q-1 1 



C(t) 



q VzTtDt 



t — ► oo 



We have introduced a new method for the analysis of 
dynamic lattice models in one dimension. This approach 
brings to light the connection between diffusion-limited 
coalescence, A + A — ► A, diffusion-limited annihilation, 
A + A — ► 0, and the g-state Potts model with Glauber 
dynamics in the zero-temperature limit. All three mod- 
els are described by the same diffusion equation; the re- 
action between particles (or kinks) is manifested in the 
boundary condition (4a), again, shared by all three mod- 
els; and the only differences arise from the different initial 
conditions implied in each case. 

For random initial conditions we obtain results which 
are in agreement with what is known from studies of 
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